Skewness as a test of the equivalence principle 
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The skewness of the large scale distribution of matter has long been known to be a probe of grav- 
itational clustering. Here we show that the skewness is also a probe of violation of the equivalence 
principle between dark matter and baryons. The predicted level of violation can be tested with the 
forecoming data from the Sloan Digital Sky Survey. 



The normalized third order moment of the galaxy dis- 
tribution, or skewness, is defined as 



S, 



< S(x) 3 > 
(< 5{x) 2 >f 



(1) 



where 8{x) is the density contrast at the point x. Its 
value at large (weakly non-linear) scales can be calculated 
exactly assuming that structure forms only via gravita- 
tional instability: in a flat universe dominated by matter 
with Gaussian initial conditions the well-known result 
is 01 S3 = 34/7. If the density contrast is smoothed 
through a window function of typical size R the skew- 
ness becomes [2, L| 



S3 = S 3 + 



d\oga 2 (R) 
dlogR 



(2) 



where a 2 is the variance of the density field smoothed 
through the same window function. In a series of pa- 
pers it has been shown that S3 remains extremely close 
to 34/7 in dark energy models Lj- LjL Ljl - in curved spaces 
0, @|, in brane- induce gravity and in Brans-Dicke 
models jjj, with deviations that hardly exceed one per 
cent in the observationally acceptable range of cosmolog- 
ical parameters. These results have shown that S3 can be 
considered one of the best probe of the gravitational in- 
stability picture at large scales [ijj- Only scenarios with 
radically different features predict values of S3 that devi- 
ate sensitively from the standard results: non-Gaussian 
initial conditions ^3], cosmic strings Cardassian cos- 
mologies [HI H3. and modified gravity models based on 
Birkhoff's lawpj. 

Since the skewness is such a good test of gravity, it 
seems interesting to ask whether it is also a good test 
of the universality of gravity, that is of the equivalence 
principle. In this paper we focus on the class of possible 
violations of the equivalence principle in which the vio- 
lating force is mediated by a scalar field. In other words, 
we investigate the effects on S3 of a scalar field coupled to 
matter in a species-dependent way. These models, first 
studied in [16], have been revived in the context of cou- 
pled dark energy (CDE), in which the same scalar field 
cou ples to matter and drives the accelerated expansion 



Let us recapitulate the calculation of the skewness as 
detailed in Q and recently reviewed in ^jj- We start by 
writing down the Newtonian equations for a pressureless 
fluid with density p, density contrast 5 and dimensionless 
peculiar velocity Vi — v pe c,i/'H , where H — aH is the 
conformal time Hubble function and a the scale factor. 
Defining the gravitational potential $ and the auxiliary 
variable A = $/(47rpa 2 ) the Newtonian equations are 



S' + Vi(l + S)vi = 0, 



v 'i + { l + ^7) v i + v j^i v i = 



-OnViA . 



(3) 

(4) 



and the Poisson equation is V^V^A = 6, where V, 
derivates with respect to comoving coordinates and the 
prime denotes derivation with respect to a = log a. 
These equations have to be complemented by the Fried- 
mann equation for W and by the matter conservation 
equations. 

We generalize now the equations introducing a scalar 
coupling to dark matter. Such a coupling is realized in 
any theory which admits in the Lagrangian a Brans-Dicke 
term of the form f(<j>,R); the low-energy limit of super- 
string theory is the most interesting example 0|. Upon 
a conformal transformation, this theory can be written 
as Einsteinian gravity in which matter and scalar field 
interact through an exchange term in their conservation 
equations (see e.g. 0,113,01): 



(c)i/;/i 



(5) 



where k 2 = 8irG and the dimensionless coupling /3(<fi) de- 
pends on the function f(<f>,R). The coupling introduces 
two distinct effects on the Newtonian equations: first, 
due to the interaction an additional force appears as a 
source; second, the matter energy density is no longer 
conserved. The first effect implies that matter feels an 
extra force due to its interaction with the scalar field 
that will add to the right-hand-side in Eq. (@J. We 
can write this term in all generality as — |f2 m r/(a, r)ViA 
where 77(0, r) is a function that in general will depend on 
time and distance. If the scalar potential is V(4>) then its 
effective mass is m 2 = d 2 V/d(j) 2 and the interaction scale 
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is A = 1/m, Then the new force introduces a Yukawa 
correction in the gravitational potential which becomes 
r](r)/r, where rj(r) = 1 + (4/3)/3 2 e~ r / A . Here, however, 
we will assume for simplicity that A is infinite, or at least 
much larger than the observed scales: on one hand, a 
very small A would be unobservable, since the interaction 
would be effectively damped; on the other, if we interpret 
the scalar field as dark energy then its interaction scale 
A is expected to be of the order of the Hubble size (see 
[2llp. It is possible however to generalize the calculations 
to a finite A . 

The second effect arises because of the non- 
conservation of the matter energy density. Eq. im- 
plies in fact in a homogeneous and isotropic metric an 
equation of the form 

p' + 3 P =-^2/3K 2 f3( ( i>)<p'p, (6) 

whose solution p = n^rriQa^e V 2 /^ J WW can ^ e m _ 
terpreted as a varying dark matter mass, p — noa~ 3 m((f>) 
with m((j)) — TriQe J 0{<i>)d<t> anc j Uq j g numeri- 

cal density of particles at present. This time-dependent 
mass introduces an extra friction term in the Euler equa- 
tions. Therefore, the scalar-Newtonian Euler equation 
for a coupled fluid can be written in the form 

< + \Fc{a)v c + (v c • ^)v c - -'^S c (a)^A c , (7) 

where the two functions, the friction F c (a) and the source 
S c (a), are in general time-dependent. In [2lJ we have 
shown that the full relativistic perturbation treatment of 
CDE reduces to Eqs. II3I7J1 in the Newtonian limit, with 

F c (a) = 2[1 + ^ - 2/jij=] , 

S c (a) = + ^ , 

where the (3 terms in F c , S c quantify the two effects due 
to the scalar interaction. Let us stress again that the 
function (3 is in general field-dependent. 

The upper bounds on a scalar interaction with baryons 
(subscript b) are very strong, of the order of /3& < 0.01 
[22]: in the following we assume that the interaction to 
baryons is effectivel y ze ro but will generalize to (3b ^= 
at the end ( in Ref. 23] it has been proposed a model in 
which such constraints can be escaped but only for suit- 
ably chosen potentials). The bounds on a coupling to 
dark matter (subscript c) are however much weaker. In 
[24| astrophysical observations were em ploy ed to derive 
(3 C < 1.5 roughly; TV-body simulations [2|| have shown 
that the dark matter halo profile depends sensitively on 
(3 C in a class of dark energy models but, due to the con- 
troversial status of the halo profile observations, it is diffi- 
cult to derive firm upper limits. Finally, in [2H we found 
that CMB requires (3 C < 0.13; however, this result as- 
sumes that the coupling remains constant throughout the 



universe lifetime and it is actually most sensitive to the 
value of (3 C at early times. Moreover, the limits obvi- 
ously depends on the assumed priors on the cosmological 
parameters, especially on the Hubble constant. We can 
summarize the observational situation with respect to (3 C 
by saying that there are no strong upper bounds to the 
present value of a scalar field coupling to dark matter; if 
(3 C varies with time then even a value of order unity at 
present is not definitively excluded. We assume therefore 
the coupling to dark matter as a free parameter and drop 
the subscript c in (3 C - Since the baryons are practically 
uncoupled, the scalar interaction violates the equivalence 
principle. The value of (3 is therefore also a measure of 
the equivalence principle violation. 

For an uncoupled and subdominant (i.e. fib <C f2 c ) 
component like the baryons Eq. becomes 

v b + ^F b {a)w b + (v h • ^)v b = -|s 6 (a)^A c , (8) 

where F b = 2 + 2H'/H, S& = fi c + fi 6 w (l c . In the stan- 
dard pure matter case Ti'/H = — 1/2 and F Cj b = S Cl b = 1, 
while in a flat universe with a scalar field component with 
equation of state — P4,/ p<p 

^ = -i[l + 3^(l-Q m )] (9) 

(here and in the following Q m = f2/j + ^b)- 

Following the notation of refs. j2] and we expand 
the scalar-Newtonian equations in a perturbation series, 
5 = J2i ^ > an d A = J2i A^. It is convenient to define 
for each component the growth function Di(a), 

5^ = D 1 (a)5^ 1) , (10) 

where 5^ is the density contrast at the initial time (as- 
sumed Gaussian distributed) and the growth exponent 
m{a) = Di'/D\. At first order we derive therefore the 
equations 

si 1) " + ^'- 3 -S b SP = 0, (11) 

4 1) " + y^'-^ 1) = °' ( 12 ) 

Asymptotically, the dark matter drives the evolution of 
the baryons, so that the two components grow with the 
same exponent m(a) but with a biased amplitude, b = 

5fr/6c. Subtracting the two equations we see that b = 
S b /[S C + m(F b - F c )} 0. In the limit of small <j>' , which 
is a typical occurrence for a slowly-rolling dark energy 
field, F c = Fb and the (anti-)bias is simply 

b= (l + 4/3 2 /3) _1 , (13) 

constant in time and space. In the same limit it appears 
that both the background evolution and the perturbation 
equations depend on (3 2 so that the sign of (3 is irrelevant. 
In it remarkable that in the opposite limit in which the 
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field kinetic energy is much larger than the potential en- 
ergy (for instance in the original Brans-Dicke model in 
which V = 0) so that the field does not drive the acceler- 
ation, it turns out that <j)' oc j3 H^] and the product (3$ 
in F c is proportional to 1 . Then even in this case the 
sign of (3 does not matter. In the following we put /3 > 0. 

Here and below we will assume for the numerical inte- 
grations an inverse power law potential V ~ 4>~ n . The 
potential appears only in the background equations and, 
indirectly, in the assumption A — > oo. For this poten- 
tial the present equation of state is approximated by 
w^o = — 2/(n + 2) |22| during the tracking regime (which 
may or may not extend to the present epoch; in the latter 
case w<po — > — 1). Integrating numerically Eqs. ll 1111211 
we find a fit for m 

m « ^ (w) 

almost independent of n (we explored the range n € 
(0,2)). 

We proceed now to second order. We define for 
each component b, c the second-order Fourier amplitude 
^ c (fc,tt) = D2b, c {ot)5\ c , {k) and, following the standard 
technique of Fourier convolution (see e.g. 0, IE El) we 
obtain for D 2b<c the equations 

D 2b + Y D *b - \s b D 2c = (^S b b + ^m 2 fc 2 ) D\{lh) 

D'L + y D 2c - \s c D 2c = (ls c + ^m 2 ) D\. (16) 

with the initial conditions D 2byC (ai n ) = D' 2b c {a in ) 
I) It is interesting to note that equations similar to 
Illlill2t have been derived in (1(1 ll5T | for a general single- 
component density expansion of Friedmann equations. 
However, the mass non-conservation induced by scalar 
gravity introduces non-standard friction terms F b , c that 
cannot be accounted for within the class explored in Ref. 

nana. nr . 

As it has been shown in |3, Ifj > the dominant term in 
the skewness is 

To derive an approximate analytical solution we can as- 
sume D 2b = b^D 2c with a constant second-order bias 
b^ since, as in the linear equations, the baryon evolution 
is driven by the dark matter one. The bias M 2 ) is not 
to be confused with the bias b 2 employed in literature 
in the Taylor expansion of a non-linear mapping from 
the underlying matter density to the galaxy distribution, 
see e.g. [IjJ. It turns out that for small (3 2 the leading 
non-trivial term is 

S 3b = & (2) 2 ( 34» m \ 

S 3c b 2 ~ P ^28m 2 + 57fi m / ' 1 ' 

(here we employed the approximation S^c w 34/7, see 
below). Substituting 111411 it appears that S 3b is almost 



independent of the potential slope n and, since m 2 « £l m , 
also almost independent of f2 m . 

The system II11I12I15I16|I . along with the background 
equations © and © , constitute a complete set of differ- 
ential equations for the unknowns 6, m, S3;,, Ss c as func- 
tion of the cosmological parameters w$,Cl m ,/3. Each of 
the functions b, m, S 3b , S 3c depends on (3 and is therefore 
in principle a test of the equivalence principle and, more 
in general, of coupled dark energy. In the long-term this 
redundancy can be exploited to set more stringent lim- 
its to the coupling and to break degeneracies with other 
cosmological parameters. However, b and S 3c require the 
detection of the large scale clustering of the dark matter 
component, while m, the clustering growth rate, requires 
accurate observations over an extended range of redshifts 
and, consequently, the problematic removal of redshift- 
dependent selection effects. Moreover, Eq. 11411 implies 
a strong level of degeneracy between (3 and the param- 
eters that enter fl m (z). On the other hand, S 3b is an 
efficient probe of the scalar interaction since it requires 
only observations of the baryon distribution at a fixed 
redshift. 

We plot in Fig. 1 the functions 53& iC (fi m ,/3) obtained 
through numerical integration. As anticipated, we find 
that S^c is close to the standard value 34/7 in the whole 
parameter range while Sa b deviates from it by more than 
1% for (3 > 0.1, following the approximate fit 

S 3b = y (1 + 0.6/? 2 )f7 m - 0005 ^ 2 ) 025 . (18) 

This result is almost independent of fl m (in fact the last 
factor can be omitted) and n and also independent of 
time (if (3 is constant) and of scale: it shows therefore 
that S3b is a direct test of the equivalence principle. 

The analytical behavior l|17H is relatively accurate (er- 
ror on 53h< 1%) only for j3 < 0.2. So far we assumed 
(3 b = but it is not difficult to see that, in the limit 
fl b -C fi C! Eqs. II14I17I18I) generalize to a finite (3 b by sim- 
ply replacing 1 with (3 C (J3 C — (3 b ). Let us remark also that 
although we performed the numerical integration with a 
dark energy potential, the scalar field need not be the 
field responsible of the accelerated expansion. The only 
condition on the potential for as concern the validity of 
our numerical calculations is that the interaction scale A 
be much larger than the astrophysical scale at which the 
observations are carried out. For instance, the original 
Brans-Dicke model, which does not give acceleration and 
where 1^ = 0, fulfills this condition (we refer here to a 
Brans-Dicke model with a species-dependent coupling as 
inRrfJIfJ). 

In [ll| the authors compiled an extensive list of present 
constraints on the smoothed skewness 5*3 from angular 
and redshift galaxy catalogs. Although several exper- 
iments quote values of S3 with errors of 5-10%, many 
results are clearly not compatible with each other. Same 
scatter, if not larger, can be seen in angular catalogs. 
This clearly points to the presence of systematic errors, 
that are likely to reside in sampling and finite volume 
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Figure 1: Contour plot of the observable quantities 536 • (7/34), 53c • (7/34) calculated numerically as function of fim, /3. For 536 
the lines correspond to the contour values 1., 1.002, 1.01, 1.03, 1.05, 1.08, while for 5 3c they are 1, 1.001, 0.999, 0.995, 0.99, 0.98, 
both from bottom to top. The short-dashed curves in the plot 536 are the fit l|18jl : the dotted curves the approximation 11711 . 
The light rectangle marks the astrophysical bounds on Sl m , the darker one adds the CMB constraint assuming constant /3 . 



effects and redshift distortions, so it is premature to per- 
form a direct comparison with data. However, analyses 
from larger redshift surveys like the Sloan Digital Sky 
Survey (SDSS) promise to measure S3 at large scales with 
a precision of less than 10% and perhaps down to 1% (see 
e.g. preliminary results in [23,[2<j] and, for the 2dF sur- 
vey, in t 30j|). At this level, SDSS might detect the scalar 
interaction or put a stringent upper limit to its present 
value. 

It is however to be stressed that our calculations refer 
to the properties of baryons, while observations deal with 
light, i.e. with the fraction of baryons that collapsed in 
sufficiently bright galaxies. The relation between the two 
populations is not well known, although at large scales, 
where hydrodynamical effects and strong non-linearities 
are smeared out, one does not expect significant segre- 
gation. To ascertain this relation it will be necessary to 
perform TV-body simulations with broken equivalence, as 
in [25] . On the other hand, it is also possible to study ob- 
jects that seem to trace with more accuracy (or at least in 



a simpler way) the underlying baryon component, such 
as Lyman-a clouds j^. Errors less than 10% in the 
bispectrum at large scales are predicted in [3]} using a 
Lyman-a forest that simulate SDSS data. 

Although models with non-Gaussian initial conditions, 
non-gravitational effects or non-standard Friedman equa- 
tions predict S3 f 34/7 [HI El El, 

they also predict a 
specific time and/or scale dependence that make them 
distinguishable, at least in principle, from a scalar inter- 
action. Further information and can be gained by the full 
bispectrum £?(ki,k2) =< 5k, ^ko^ -k, -k, >, rather than 
by the integrated skewness. In [23] it has been shown that 
the scale-dependence of the bispectrum may be of great 
help in constraining primordial non-gaussianity. The be- 
havior of the bispectrum for the present model will be re- 
ported in subsequent work. Forecoming large scale skew- 
ness data offer therefore the exciting opportunity to test 
the equivalence principle in a realm inaccessible to labo- 
ratory experiments. 
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